Chapter 3 



Collisions in Plasmas 



3.1 Binary collisions between charged particles 

Reduced-mass for binary collisions: 

Two particles interacting with each other have forces 

Fi2 force on 1 from 2. 

F21 force on 2 from 1. 

By Newton's 3rd law, F12 = — F21. 

Equations of motion: 



Combine to get 
which may be written 



miTi = F12 ; m2r2 = F21 (3.1) 



ri - r2 = F12 (— + — ] (3.2) 



(ri - r2) = F12 (3.3) 



mi + 777.2 dt'^ 

If F12 depends only on the difference vector ri — r2, then this equation is identical to the 
equation of a particle of "Reduced Mass" 777^ = moving at position r = ri — r2 with 

respect to a fixed center of force: 

777,r = Fi2(r) . (3.4) 

This is the equation we analyse, but actually particle 2 does move. And we need to recognize 
that when interpreting mathematics. 

If F21 and ri — r2 arc always parallel, then a general form of the trajectory can be written 
as an integral. To save time we specialize immediately to the Coulomb force 

= 4 (3.5) 

Solution of this standard (Newton's) problem: 
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1: projectile 

Vi 



b Impact Parameter' 

2: target 

Figure 3.1: Geometry of the collision orbit 



Angular momentum is conserved: 

rrirr'^O = const. = rrirbvi 

Substitute u=l then 9 = ^ = u'^ bvi 
Also 

d 1 



clockwise from symmetry) 



1 du 



dt u dO 
d^u ■ 



-bvi 



du 

de 



-bvi 



d9^ 



dO^ 



Then radial acceleration is 



f — rO^ = — (6f i)^ 



+ u 



'12 



I.e. 



d'^u 



+ u 



Q1Q2 



dO'^ ' 47reo rrir (bvif 

This orbit equation has the elementary solution 



u = - = C cos u — 



1 



47reo nir {bviY 



(3.6) 



(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 



The sin^ term is absent by symmetry. The other constant of integration, C, must be deter- 
mined by initial condition. At initial (far distant) angle, 6*1, mi = — = 0. So 



= Ccos^i 



47reo nir {bviY 



There: 



Hence 



Ti = —Vi = —bvi — \i = +bviC sm&i 
dO 



tan 6*1 



-1/Cb 



sin 61 

cos^i " ,1 ,J C 

^ Atteo mr{bvi) ' 



b 

bgo 



(3.12) 

(3.13) 
(3.14) 
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where 

090 = 2 ■ (3-15) 

47reo nirVf 

Notice that tan6'i = — 1 when b = 699. This is when 9i = —45° and x = 90°- So particle 
emerges at 90° to initial direction when 

b = 690 "90° impact parameter" (3.16) 

Finally: 

C^--cosec9, = --[l + f'j (3.17) 
3.1.1 Frames of Reference 

Key quantity we want is the scattering angle but we need to be careful about reference 
frames. 

Most "natural" frame of ref is "Center-of-Mass" frame, in which C of M is stationary. C of 
M has position: 

mi + m2 



and velocity (in lab frame) 



Now 



mi + m2 



ri = R+ r (3.20) 

mi + 777.2 

r2 = R — — r (3.21) 

mi + 7772 

So motion of either particle in C of M frame is a factor times difference vector, r. 
Velocity in lab frame is obtained by adding V to the C of M velocity, e.g. ^^^^ + V. 
Angles of position vectors and velocity differences are same in all frames. 
Angles (i.e. directions) of velocities are not same. 

3.1.2 Scattering Angle 

In C of M frame is just the final angle of r. 

- 2^1 + X = TT (3.22) 

(^1 is negative) 

X^7^ + 20i ; ^1 = ^^ ■ (3.23) 
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Figure 3.2: Relation between 6i and %. 



tan 9i = tan f^ = ^\=- cot | (3.24) 

So 

cotf = ^ (3.25) 

2 090 

tan| = ^ (3.26) 

But scattering angle (defined as exit velocity angle relative to initial velocity) in lab frame 

is different. 

Final velocity in CM frame 

v'cM = T;icM(cosxc,sinxc) = — — t-i (cos Xc, sin Xc) (3.27) 

mi + m2 

[ Xc = X and v\ is initial relative velocity] . Final velocity in Lab frame 

= vU + V = + cos xe, sm xe) (3.28) 

V mi + m2 mi + m2 / 

So angle is given by 

cot XL = ,„ ' . = cosecxc + cot Xc 3.29) 

, ^ sm y,,. v^ mi 

For the specific case when m2 is initially a stationary target in lab frame, then 

V — ^^^''^ — and hence (3.30) 

mi + m2 

cotxL = — ^cosecXc + cot Xc (3.31) 
m2 

This is exact. 

Small angle approximation (cotx ^^cosecx — ^ gives 

1 mi 1 1 m2 

— = —- + - ^ XL = 3.32 

XL rn2 Xc Xc mi + m2 

So small angles are proportional, with ratio set by the mass-ratio of particles. 
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Figure 3.3: Collisions viewed in Center of Mass and Laboratory frame. 



3.2 Differential Cross-Section for Scattering by Angle 

Rutherford Cross-Section 

By definition the cross-section, a, for any specified collision process when a particle is passing 
through a density n2 of targets is such that the number of such collisions per unit path length 
is n2(T. 

Sometimes a continuum of types of collision is considered, e.g. we consider collisions at 
different angles (x) to be distinct. In that case we usually discuss differential cross-sections 
(e.g ^) defined such that number of collisions in an (angle) element dx per unit path length 
is n2^dx. [Note that ^ is just notation for a number. Some authors just write cr(x), but I 
find that less clear.] 

Normally, for scattering-angle discrimination we discuss the differential cross-section per unit 
solid angle: 

da 

P.33) 

This is related to scattering angle integrated over all azimuthal directions of scattering by: 
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Figure 3.4: Scattering angle and solid angle relationship. 



dQs = 27rsinxdx (3.34) 

So that since 

da ,^ da , , , 

dQs = -r'^X (3.35) 



dVts dx 
we have 

da \ da , , 

dVts 2-7? sin X dx 

Now, since x is a function (only) of the impact parameter, 6, we just have to determine the 
number of collisions per unit length at impact parameter h. 




Figure 3.5: Annular volume corresponding to dh. 



Think of the projectile as dragging along an annulus of radius h and thickness dh for an 
elementary distance along its path, dl. It thereby drags through a volume: 

dl2Tihdh . (3.37) 

Therefore in this distance it has encountered a total number of targets 

d^2^lhdh . n2 (3.38) 
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at impact parameter b(db). By definition this is equal to M^dbn2. Hence the differential 
cross-section for scattering (encounter) at impact parameter b is 

27r6 . 



da 
db 



Again by definition, since x is a function of b 



da 
dx 



dx 



da 
db' 



-db 



da 
dx 



da 
db 



db 
dx 



[db/dx is negative but differential cross-sections are positive.] 
Substitute and we get 



da 
dQ^ 



1 da 
27r sin x db 



db 
dx 



b 

sin 



db 
dx 



[This is a general result for classical collisions.] 
For Coulomb collisions, in C of M frame. 



db 
dx 



= b 



'90 



cot 

d 

dx 



b 

ho 



. X 
cot — — 
2 



090 2% 

— cosec — 
2 2 



Hence 



da 
dQg 



'90 cot f 690 



2% 

— cosec — 
2 2 



smx 

2 

-'90 

2 2 sin I cos | sin^ | 



foon cos f / sin I 



bio 



4sin4| 



This is the Rutherford Cross-Section. 



da 
dfl. 



b"^ 



4sin^f 



(3.39) 



(3.40) 



(3.41) 



(3.42) 
(3.43) 

(3.44) 
(3.45) 
(3.46) 



(3.47) 



for scattering by Coulomb forces through an angle x measured in C of M frame. 
Notice that — > 00 as x ^ 0. 

This is because of the long-range nature of the Coulomb force. Distant collisions tend to 
dominate, (x ^ <^ 6 ^ 00) . 
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3.3 Relaxation Processes 



There are 2 (main) different types of collisional relaxation process we need to discuss for a 
test particle moving through a background of scatterers: 

1. Energy Loss (or equilibrium) 

2. Momentum Loss (or angular scattering) 

The distinction may be illustrated by a large angle (90°) scatter from a heavy (stationary) 
target. 

If the target is fixed, no energy is transferred to it. So the energy loss is zero (or small if 
scatterer is just 'heavy'). However, the momentum in the x direction is completely 'lost' in 
this 90° scatter. 

This shows that the timescales for Energy loss and momentum loss may be very different. 
3.3.1 Energy Loss 

For an initially stationary target, the final velocity in lab frame of the projectile is 

, f mivi m2Vi m2Vi . \ 

Vt = \ cosxc , smxc (3.48) 

\mi + 777.2 n^i + ^^2 iT^i + ' 

So the final kinetic energy is 

1 12 1 2// 2mi7772 



= ^miv'l = ^mivt\ [ +7 cosXc (3.49) 

2 2 [Vr?7i + 7772/ (mi + 7772) 



+ 7 ""^^ (cos'xc + sin^xc) \ (3.50) 

(mi + m2) ^ ^ J 

= ^mi^;i^ 1 1 + ^"^^^^ (cosXc-l)[ (3.51) 
^ t (mi + m2) J 



-mi^i < i + ^ 

^ I (mi + m2 



mi^i < 1 H 2 2sm — > (3.52) 



Hence the kinetic energy lost is AK — K — K' 

1 4mim2 . 2 Xc .0 r,o\ 

= ^m^V2 -2 sm — (3.53) 

^ (mi + m2j ^ 

1 2 4mim2 1 . . ,Xc h 

— -miV^ -2 2 [usmg cot — = —J (3.54) 

2 (mi + m2) ( A-\ +1 2 O90 

(exact). For small angles x ^ 1 i-e- b/bgo ^ 1 this energy lost in a single collision is 
approximately 



1 2\ 4mim2 (bgo 



2 



mi.n -^=2 X (3.55) 
^ / (mi + m2) V " / 
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If what we are asking is: how fast does the projectile lose energy? Then we need add up the 
effects of all collisions in an elemental length M at all relevant impact parameters. 

The contribution from impact parameter range dh at h will equal the number of targets 
encountered times A.K: 

^eo h^h ^ 2 4mim2 (bgoV 
n2di27ibab -niiv^ 2 — (3.56) 



encounters 



2 [nil + 1112) \ h 



Loss per encounter (AK) 
This must be integrated over all b to get total energy loss. 

l/;/i//i2 f fbgo^'^ 



dK = n2dlK '—^ \^\ ^T^bdb (3.57) 

(mi + 7712) \ b J 

so 

§ = ^ -2-^^^ 8nbl, [In \b\]Z: (3.58) 
«t (mi + m2) 

We see there is a problem both limits of the integral [b ^ 0, b ^ 00) diverge logarithmically. 
That is because the formulas we are integrating are approximate. 

1. We are using small-angle approx for AK. 

2. We are assuming the Coulomb force applies but this is a plasma so there is screening. 

3.3.2 Cut-offs Estimates 

1. Small-angle approx breaks down around b — bgo- Just truncate the integral there; 
ignore contributions from b < bgo- 

2. Debye Shielding says really the potential varies as 



cxp [t- 1 

(j) oc instead of oc - (3.59) 

r r 

so approximate this by cutting off integral at 6 = equivalent to 

brain = ^max = Aq. (3.60) 

f = A-„.-^^^8.6^1n|A| (3.61) 

at (mi -|- m2) 

Ad feoTeY / ( _QiQ2_ 



^ ^ (3_g2) 

690 \ne^ J / \4:7TeomrVfl 
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So Coulomb Logarithm is 'In A' 

A 



epTe 

2 



91 92 



Because these cut-ofFs are in In term result is not sensitive to their exact values. 
One commonly uses Collision Frequency. Energy Loss Collision Frequency is 



Substitute for 690 and (in 690) 



1 dK mim2 „ ,2 , 1.1 

navi -2 877690 In |A| 

(mi + 777,2) 



n2Vi 



mim2 



Stt 



Stt 



91 92 



' m\+m2 



In A 



7l2 



9i9i 



(47reo)^ mim2vf 
Collision time = I/i^k 

Effective (Energy Loss) Cross-section 



In A 



1 dK 
K <U 



0-^712 



Stt 



(47reo) 77717772^1 



In A 



(3.63) 



(3.64) 



(3.65) 
(3.66) 



(3.67) 



3.3.3 Momentum Loss 

Loss of x-momentum in 1 collision is 



= 777,1(^1 

= 777,i7;i |l - 
7772 

= Px 



-'Lxj 



mi 



7772 



777i + 7772 



+ 

777i + 7/72 7771 + 7772 
(1 -COSXc) 



COSXc 



)} 



Px 



7772 



777i + 7772 2 



Px 



7772 



2^2 



777i + 7772 



(3.68) 

(3.69) 
(3.70) 

(3.71) 



(small angle approx). Hence rate of momentum loss can be obtained using an integral 
identical to the energy loss but with the above parameters: 



dp 
d£ 



7l2P 
772P 



m2 



bmin 7771 + 7772 
^2 , ,2 



2b^ 

2Mb 



mi + 7772 



47r In A 



Note for the future reference: 

dp 



dp 2 777 i7r72 

dt di mi + 7772 



inblg In A. 



(3.72) 
(3.73) 

(3.74) 
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Therefore Momentum Loss. 
Collision Frequency 

1 dp 1712 



Vr, = vi-—^n2Vi — — A-Kh^Q In A (3.75) 

n 2 



p (M. nil + 



7112 . 

n2Vi 47r 

mi + m.2 



"mi+jTt2 -I 



In A (3.76) 



^ n ^1^2 47r (mi + ma) 
^ (47reo)^ m2mlvl 



Collision Time Tp = l/i^p 
Cross-Section (effective) a = i^p/n2Vi 
Notice ratio 



Energy Loss z/k 2 /mi + m2 2mi 



Momentum loss i/p mima / mamf mi + m2 

Tfiis is 



(3.78) 



~ 2 if mi » m2 (3.79) 
= 1 if mi = ma (3.80) 
< 1 if mi < m2. (3.81) 

Third case, e.g. electrons shows that mostly the angle of velocity scatters. Therefore 
Momentum 'Scattering' time is often called '90° scattering' time to 'diffuse' through 90° in 
angle. 

3.3.4 'Random Walk' in angle 

When mi << m2 energy loss << momentum loss. Hence |v^ ~ t^i. All that matters is the 
scattering angle: xl — Xc — '^bgo/b. 

Mean angle of deviation in length L is zero because all directions are equally likely. 

But: 

Mean square angle is 



Aq;2 = n2L / X '^T^hdh (3.82) 
= Ln2 87r hl^ In A (3.83) 



Spread is 'all round' when Aq;^ ~ 1. This is roughly when a particle has scattered 90° on 
average. It requires 

Ln2 87r6^o lnA = 1 . (3.84) 
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So can think of a kind of 'cross-section' for 'ergo' 90° scattering as such that 

n2l/'cr9o' = 1 when Ln2 Stt 6go InA = 1 (3.85) 
i.e. 'cTgo' = STrftgolnA (= 2(7p) (3.86) 
This is 8 In A larger than cross-section for 90° scattering in single collision. 

Be Careful! 'cxgo' is not a usual type of cross-section because the whole process is really 
diffusive in angle. 

Actually all coUision processes due to coulomb force are best treated (in a Mathematical 
way) as a diffusion in velocity space 

— > Fokker-Planck equation. 
3.3.5 Summary of different types of collision 

The Energy Loss collision frequency is to do with slowing down to rest and exchanging 
energy. It is required for calculating 

Equilibration Times (of Temperatures) 
Energy Transfer between species. 
The Momentum Loss frequency is to do with loss of directed velocity. It is required for 
calculating 

Mobihty : Conductivity / Resistivity 

Viscosity 

Particle Diffusion 

Energy (Thermal) Diffusion 
Usually we distinguish between electrons and ions because of their very different mass: 

Energy Loss [Stationary Targets] Momentum Loss 



K 



K 



K 



Vi, 



57r 

\2 ,™2«i3 



u^i = n.. 



(47reo) rnlv^ 
Z'^e^ 87r 
(47reo)^ meruiv, 



In A 



In A "u,. 



K 



In A 



me + me 
. 2me 
me + rui 



rrij 
2mp 



1 



TT 



In A Pvi^^^Vi, 



X 



2mi 
me + rrii 1 
2mi- 2 



(3.87) 



Sometimes one distinguishes between 'transverse diffusion' of velocity and 'momentum loss'. 
The ratio of these two is 





dxl 1 


1 d-p 




dL 1 


p dL 




( mi 






Vmi+m2 









mi+m2 2 



2m2 
mi + m2' 



(3.88) 
(3.89) 
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So 



Hence 



'cTgo' 2m2 



'cTp' mi + 7712 



= 1 like particles (3.90) 

~ 2 mi « m2 (3.91) 

~ — m2 « mi . (3.92) 
mi 



= Ve= Ve(= Vee'!!) (3.93) 

^l/ei = 2fl/ei= Ve-^'(=^i^ee)(= Vei') (3.94) 

Tig 

-^i/jj = Pj^ii = ^Uii (= z/jj!!) (Like Ions) (3.95) 

!_ 2me „ rrif. k k /q 

[But note: ions are slowed down by electrons long before being angle scattered.] 



3.4 Thermal Distribution Collisions 

So far we have calculated collision frequencies with stationary targets and single-velocity 
projectiles but generally we shall care about thermal (Maxwellian) distributions (or nearly 
thermal) of both species. This is harder to calculate and we shall resort to some heuristic 
calculations. 



3.4.1 



Very rare for thermal ion velocity to be ~ electron. So ignore ion motion. 
Average over electron distribution. 

Momentum loss to ions from (assumed) drifting MaxweUian electron distribution: 

\2- 



2T 



(3.97) 



Each electron in this distribution is losing momentum to the ions at a rate given by the 
collision frequency 

q^Qi Att (me + rrii 



rii 



(47reo)^ 



m^mg^;^ 



In A 



so total rate of loss of momentum is given by (per unit volume) 

dp 
~ dt 

To evaluate this integral approximately we adopt the following simplifications. 



j /e(v) iyp{v) nieV d^v 



(3.98) 



(3.99) 
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1. Ignore variations of In A with v and just replace a typical thermal value in A = 

><D/b9o{vi). 

2. Suppose that drift velocity is small relative to the typical thermal velocity, written 
Ve = ^Jtjrrie and express /g in terms of u = ^ to first order in = 



fe 



Tie 



(27r) 2 V, 



exp 



-1 f 



Up 



[1 + u.Ud) exp 



-u 



[1 + U:cUd) fo 



taking x-axis along and denoting by fo the unshifted Maxwellian. 
Then momentum loss rate per unit volume 



dPx [ f ,3 



dt 



/V 
{l + UxUd)fo ^ Vx V 



To evaluate this integral, use the spherical symmetry of fo to see that: 



6 Jo V 

rex 

/ /oS 

Jo 



3 
27r 



r2vdv 



3 (27r)5i;3^o 
2n He 2 _ _ 
3 (27r)5 



exp 



— f 
2^ 



r "-e 

2 



Thus the Maxwell-averaged momentum-loss frequency is 

1 dp 2 
p dt 



3(27r 



3(27r 

(where p = meVdUe is the momentum per unit volume attributable to drift). 



r 

3(27r)2 (47reo 



miimlvl 



In A, 



2 / ze^ \ ^ 47r 



3(27r) 



,47ren 



InAg 



(3.100) 
(3.101) 



(3.102) 



(3.103) 



(3.104) 



(3.105) 



(3.106) 



(substituting for thermal electron velocity, I'g, and dropping ^ order term), where Ze — qi. 
This is the standard form of electron collision frequency. 
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3.4.2 i e 



Ion momentum loss to electrons can be treated by a simple Galilean transformation of the 
e ^ i case because it is still the electron thermal motions that matter. 

Ions — > Electrons ions < — Electrons 




Vi v^ = ^ Vd = -Vi 



Figure 3.6: Ion-electron collisions are equivalent to electron-ion collisions in a moving refer- 
ence frame. 



Rate of momentum transfer, is same in both cases: 

-pu (3.107) 



dp 
dt 

Hence pei^ei = Pi^^ie or 



Vie = —l^ei = J^ei (3.108) 

Pi niUli 

(since drift velocities are the same). 

Ion momentum loss to electrons is much lower collision frequency than e ^ i because ions 
possess so much more momentum for the same velocity. 



3.4.3 i ^ i 

Ion-ion collisions can be treated somewhat like e ^ i collisions except that we have to 
account for moving targets i.e. their thermal motion. 

Consider two different ion species moving relative to each other with drift velocity v^, the 
targets' thermal motion affects the momentum transfer cross-section. 

Using our previous expression for momentum transfer, we can write the average rate of 
transfer per unit volume as: [see 3.74 "note for future reference"] 

I l^r vAn bl, InA hhd\,d\, (3.109) 

dt J J mi + 1712 



where is the relative velocity (vi — V2) and 690 is expressed 



-"90 



47760 rrij-vl 
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(3.110) 



and rrir is the reduced mass '^^J^^ . 

' mi+m2 



Since everything in the integral apart from /1/2 depends only on the relative velocity, we 
proceed by transforming the velocity coordinates from Vi, V2 to being expressed in terms of 
relative (v^) and average (V say) 



Vi - V2 



V 



miVi + m2V2 



mi + 777.2 

Take /i and /2 to be shifted Maxwellians in the overall C of M frame: 

\2n 



(3.111; 



( ^3 \ 
\2t:T) 



fj = ^3 

where miVdi + m2Vd2 — 0. Then 



exp 



2T 



/1/2 



nin2 



^ rrii \ 2 / 1112 \ 
\2nT) UttT; 



exp 



(J = l,2) 



(3.112) 



2T 



2T 



(3.113) 



to first order in v^. Convert CM coordinates and find (after algebra) 

= ^^^^ (2^) ^ (^) ^ 



rrirV^ 
~2T 2T~ 



X U + ^Vrf.V, 



(3.114) 



where M — mi + m2. Note also that (it can be shown) dPvi(fv2 — (fvrdPV. Hence 



dp 



II 



VrTTlrVr^'K 690 In An 1^2 



M \-2 f m. 



2'kT) \2'kT) 



exp 



2T 



exp 



rrirvl 
2T 



^[l + '^va.Vr]d\£V (3.115) 



and since nothing except the exponential depends on V, that integral can be done: 
rrirVr'iTr In A nin2 



dp 

~dt 



f rnr_\ 
\2nT) 



exp 



27r 



rrir 



l + ^Vd.v,^dV (3.116) 



This integral is of just the same type as for e — i coUisions, i.e. 



dp 
dt 



VdVrtrrir 47r hg^ivrt) In At nin2 / —fo{vr)d-^Yr 

U3 



VdVrtrrirA TT 690 (t'rt) In At nin2 



3(27r) 



(3.117) 
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where Vrt = y ^90 ('^rt) the ninety degree impact parameter evaluated at velocity Vtr, 
and fo is the normalized Maxwellian. 

- — = r hi hiKtnin2mrVd (3.118) 

3(27r)2 V47reo/ m^w^ 

This is the general result for momentum exchange rate between two Maxwellians drifting at 
small relative velocity Vd. 

To get a collision frequency is a matter of deciding which species is stationary and so what the 
momentum density of the moving species is. Suppose we regard 2 as targets then momentum 
density is niUiiVd so 



_ 1 dp 

^12 = 



r".fl^l'^— . (3.119) 

2 \47reo/ mrV^i mi 



nimiVddt 3(27r)2 V47reo/ rrirV^ 

This expression works immediately for electron-ion collisions substituting — mg, recov- 
ering previous. 

2 / I 

For equal- mass ions rrir = — = Irrii and Vrt = \ — = \ —■ 
Substituting, we get 

T^ii = A (7^)' In A (3.120) 

37r2 \47reo/ ^22^2 

that is, ^ times the e — i expression but with ion parameters substituted. [Note, however, 
that we have considered the ion species to be different.] 



3.4.4 e ^ e 

Electron-electron collisions are covered by the same formalism, so 



1 / eM^ 47r 
V,, = -^ne[- -^-^ InA . (3.121) 

However, the physical case under discussion is not so obvious; since electrons are indistigu- 
ishable how do wc define two different "drifting maxwellian" electron populations? A more 
specific discussion would be needed to make this rigorous. 

Generally u^e ~ ^ei / \/2 : electron-electron collision frequency ~ electron-ion (for momentum 
loss) . 

3.4.5 Summary of Thermal Collision Frequencies 

For momentum loss: 

Vei = ^n,[^^ InAe . (3.122) 
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i/ge — ^ei ■ (electron parameters) (3.123) 

_^ ^ ^j,^^ _ (3.124) 
niUii 

Fjj/ = ^ rii, ( \ — j-^^ ( — — \ InAj (3.125) 

3V7r V47reo/ ^sj^s Vrnj + mj,/ 



Energy loss related to the above {^u) by 



mi + 1712 



Transverse 'diffusion' of momentum ^u, related to the above by: 



-jy^^I^Pu . (3.127) 

mi + 777.2 



3.5 Applications of Collision Analysis 
3.5.1 Energetic ('Runaway') Electrons 

Consider an energetic {^rrievl » T) electron travelling through a plasma. It is slowed down 
(loses momentum) by collisions with electrons and ions (Z), with collision frequency: 



Vee = Uee^Ue— -3 In A (3.128) 



"l^ei = ^ZUee (3.129) 



Hence (in the absence of other forces) 

d 



^^{nieV) = - {^Uee + ^I'ei) iripV (3.130) 

Z\ 

1 + — ) fee "^ef (3.131) 



2 

This is equivalent to saying that the electron experiences an effective 'Frictional' force 

d ( Z \ 

= Tr("^e^^) = - (1 + 77 ) '^ee (3.132) 



dV ^ ' V 2 

Z\ Stt In A 

— I Tl 

2; " (47reo)' m^v-^ 



Ff = -(l+o)^e 7^.37^—^ (3.133) 
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Notice 

1. for Z = 1 slowing down is | on electrons | ions 

2. Ff decreases with v increasing. 

Suppose now there is an electric field, E. The electron experiences an accelerating Force. 
Total force 

d, ^ / Z\ STTlnA , , 

F = —(mv) = -eE + Ff ^ -eE - 1 + - ^ 3.134 

Two Cases (When E is accelerating) 

1. \eE\ < \Ff\: Electron Slows Down 

2. |e£;| > |F/|: Electron Speeds Up! 

Once the electron energy exceeds a certain value its velocity increases continuously and the 
friction force becomes less and less effective. The electron is then said to ahve become a 



runaway . 
Condition: 



-ruev^ > 1 + - ne n — 3.135 

2 ' V 27 ' (47reo)' 2eE ^ ' 



3.5.2 Plasma Resistivity (DC) 

Consider a bulk distribution of electrons in an electric field. They tend to be accelerated by 
E and decelerated by coUisions. 

In this case, considering the electrons as a whole, no loss of total electron momentum by 
e — e collisions. Hence the friction force we need is just that due to Pgi- 

If the electrons have a mean drift velocity Vci{« Vthe) then 

^{rrieVd) = -eE - V^im^Vd (3.136) 
at 

Hence in steady state 

Vd=^^ . (3.137) 

n^e'^E 

j = -rieevd = — — (3.138) 
Now generally, for a conducting medium we define the conductivity, a, or resistivity, rj, by 

j = aE ; rij = E ('^ = ^) (3-139) 
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The current is then 



Therefore, for a plasma, 

a = ' = ^ (3.140) 

Substitute the value of I^ei and we get 

mZ^ e^ml Stt In A 
7] ~ • 3- (3.141) 

(47reo)' 3v^Te^ 
Ze^ml Stt In A 



(for a single ion species). (3.142) 

(47reo)' 3V27rre^ 



3 

^ e 



Notice 

1. Density cancels out because more electrons means (a) more carriers but (b) more 
collisions. 

2. Main dependence is 77 oc T~'^l'^ . High electron temperature implies low resistivity (high 
conductivity) . 

3. This expression is only approximate because the current tends to be carried by the 
more energetic electrons, which have smaller Vf.i\ thus if we had done a proper average 
over /(I'e) we expect a lower numerical value. Detailed calculations give 

?7 = 5.2 X 10"^ 3 (3.143) 

for Z = \ (vs. ~ 10~^ in our expression). This is 'Spitzer' resistivity. The detailed 
calculation value is roughly a factor of two smaller than our calculation, which is not 
a negligible correction! 



3.5.3 Diffusion 

For motion parallel to a magnetic field if we take a typical electron, with velocity v\\ ~ Vte it 
will travel a distance approximately 

4 = Vte/lye^ (3.144) 

before being pitch-angle scattered enough to have its velocity randomised. [This is an order- 
of-magnitude calculation so we ignore Uee-] ^ is the mean free path. 

Roughly speaking, any electron does a random walk along the field with step size i and step 
frequency Uei- Thus the diffusion coefficient of this process is 

Deiic^eluei^'^ . (3.145) 

Similarly for ions 

All ^ ipii ^ ^ (3.146) 
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Notice 

v^lv^ ^["^^^ ^""-^ (if ^ T,) (3.147) 

Hence 4 — 4 

Mean free paths for electrons and ions are ~ same. 
The diffusion coefficients are in the ratio 

D fm \ i 

— ^ ~ I — - j : Ions diffuse slower in parallel direction. (3.148) 
Diffusion Perpendicular to Mag. Field is different 





Figure 3.7: Cross-field diffusion by collisions causing a jump in the gyrocenter (GC) position. 

Roughly speaking, if electron direction is changed by ~ 90° the Guiding Centre moves by 
a distance ~ r/,. Hence we may think of this as a random walk with step size ~ and 
frequency z7ei- Hence 

Dex ^ rl,V,, ~ ^F,, (3.149) 

Ion transport is similar but requires a discussion of the effects of like and unlike collisions. 

Particle transport occurs only via unlike coUisions. To show this we consider in more detail 
the change in guiding center position at a coUision. Recall mv = gv A B which leads to 

v_L = — r^ A B (perp. velocity only). (3.150) 
m 

This gives 

At a collision the particle position does not change (instantaneously) but the guiding center 
position (ro) does. 

r[, + v'l = ro + ri =^ Aro = r[, - ro = -(r'^ - r^,) (3.152) 
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Change in is due to the momentum change caused by the coUision: 

r'i - rz, = ^ A m(v' - v) = ^ A A(mv) (3.153) 

So 

Aro = A A(mv). (3.154) 

The total momentum conservation means that A(mv) for the two particles colliding is equal 
and opposite. Hence, from our equation, for like particles, Aro is equal and opposite. The 
mean position of guiding centers of two colliding like particles (roi +ro2)/2 does not change. 

No net cross field particle (guiding center) shift. 

Unlike collisions (between particles of different charge q) do produce net transport of particles 
of either type. And indeed may move roi and ro2 in same direction if they have opposite 
charge. 

_ _ 

A± ^ rl, fl/ie - ^ Pl^ie (3.155) 
Notice that rljrl^ ~ rrii/me ; P^ie/i^ei - ^ 

So Di±/ De± ~ 1 (for equal temperatures). CoUisional diffusion rates of particles are same 
for ions and electrons. 

However energy transport is different because it can occur by like-like coUisins. 
Thermal Diffusivity: 

Xe ~ (pei + l^ee) ~ ^^le T^ei (pei ~ l^ee) (3.156) 

Xi ~ rlii^ie + T^ii) -rliTyu {yii»Vi^) (3.157) 



.2 T7 ^ ^2 



m,- mi I rrij \ 2 



Xi/Xe ~ Lf^^^[!^=ii::iy (equal T) (3.158) 

CoUisional Thermal transport by Ions is greater than by electrons [factor ~ (mi/me)^ ]■ 
3.5.4 Energy Equilibration 

If Tg 7^ Ti then there is an exchange of enegy between electrons and ions tending to make 
Te — Ti. As we saw earher 

K 2r?T.e p Trig , . 

Uei = ^Pei = l^ei (3.159) 



So applying this to averages. 



^Vei ^ l^ei (- J^ie) (3.160) 



rrii 
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Thermal energy exchange occurs ~ mg/mj slower than momentum exchange. (Allows Tg ^ 
Ti). So 

"^^^ "^^^ -^e.(Te-r,) (3.161) 



- e 



Prom this one can obtain the heat exchange rate (per unit volume), Hgi, say: 

= --n-in-Ti)^ -n^,i{Te-Ti) (3.163) 



Important point: 

'Electrons and Ions equihbrate among themselves much faster than with each other' 



^UeiC^'^Zu^e^^J^Yuu. (3.164) 



3.6 Some Orders of Magnitude 

1. In A is very slowly varying. Typically has value ~ 12 to 16 for laboratory plasmas. 

2. i7ei 6 X 10-11 (ni/m^) / {TjeVf^ (In A = 15, Z = 1). 

e.g. = 2 X lO^fi-i (when n = lO^^m^^ and Tg = IkeV.) For phenomena which happen 
much faster than this, i.e. r -C l/t'ei ~ S/xs, collisions can be ignored. 
Examples: Electromagnetic Waves with high frequency. 

3. Resistivity. Because most of the energy of a current carrying plasma is in the B field 
not the K.E. of electrons. Resistive decay of current can be much slower than I/gj. E.g. 
Coaxial Plasma: (Unit length) 

Inductance L — jio In ^ 
Resistance R — rj I/tto^ 
L/R decay time 

yUoVra^ b n^e"^ ^ ^ 

tr ^ m - ~ — ^/^o ^ct In ~ 

77 a rripi'f.i a 



n^e^ a? 1 cj^a^ I 1 



» — . (3.165) 



Comparison 1 keV temperature plasma has ~ same (conductivity/) resistivity as a slab 
of copper (~ 2 X lO^^fim). 

Ohmic Heating Because rj oc T^g if we try to heat a plasma Ohmically, i.e. by simply 
passing a current through it, this works well at low temperatures but its effectiveness falls 
off rapidly at high temperature. 
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Result for most Fusion schemes it looks as if Ohmic heating does not quite yet get us to the 
required ignition temperature. We need auxiUiary heating, e.g. Neutral Beams. (These slow 
down by collisions.) 
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